Revisiting renormalization schemes in a differential equation approach 
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We reconsider the choice of renormalization schemes in a differential-equation approach to aid 
the discussion of the renormalization of the unstable particles and the CKM matrix in the Standard 
Model. Certain mass dependent schemes do not satisfy these natural differential equations modulo 
trace anomaly. By the way, the Callan-Symanzik equations were employed to show that both mass 
dependent and independent schemes realize fermion decoupling in the same way. 



Theoretically, different renormalization schemes arc perfectly equivalent provided the perturbation is completed. 
The perturbative scheme dependence is due to the truncation of the perturbation series [1]. Recently, the issue of 
renormalization scheme has been intensively investigated in Standard Model (SM) and its extensions with respect 
to the unstable sector (like W ± , Z° and Higgs particles) [2,3] and that of the Cabibbo-Kobayashi-Maskawa (CKM) 
' matrix [4,5], where certain mass dependent (like MOM and on-shell) schemes become problematic in these contexts: 
^ . the on-shell schemes (mass dependent) would bring in severe gauge dependence and IR singularities or ambiguities. 
The lesson that could be drawn from these progresses is that the freedom in choosing renormalization prescriptions 
deserves closer investigation in applications. In general, a gauge invariant definition is preferred while gauge dependent 
ones should be avoided at all [3,6]. 

In this brief report, we wish to draw attention to other aspects of the same issue. We start with the standard point 
of view that all the known QFT's are effective theories for the low energy sectors in a completely well-defined quantum 
, theory containing 'correct' high energy details [7]. Then we can derive from this postulate a set of natural differential 
' equations for the loop amplitudes due to the well-known fact that differentiating the loop amplitudes with respect 
to external parameters like momenta and mass(es) reduces the divergence degree 1 : the divergent loop amplitudes 
are therefore obtained as solutions to these well-defined differential equations, which are generally finite nonlocal 
terms plus certain a polynomial in terms of external momenta, masses and couplings with arbitrary coefficients to be 
fixed through physical 'boundary conditions' [8] as we usually do in quantum mechanics and electrodynamics. We 
have demonstrated this simple approach in a recent study of dynamical symmetry breaking in the two loop effective 
potential of masslessA</> 4 [9]. Obviously, for convergent diagrams such differential equations are naturally satisfied. 
The key point is that, the coefficients in the polynomial, to be fixed through physical 'boundary conditions', should 
O |' be independent of the differentiation parameters like masses. If one affords some intermediate unphysical definition 
to these constants, one could arrive at any renormalization prescription, like MOM, on-shell and MS. However, 
some prescription might be inconsistent: the original differential equations might be violated by such prescriptions. 
Since such differential equations are well defined, the physical amplitudes must satisfy these equations and any 
renormalization prescriptions that solve these equations could be used. However, it is unclear how those prescriptions 
that fail to satisfy such equations could finally lead to physical amplitudes that satisfy the same equations, as a 
practical complete summation of the perturbations series is impossible. 

Now let us demonstrate the differential equation solution of a loop amplitude with a simple vertex function at the 
lowest order in QED, say, the I-loop photon vacuum polarization tensor 
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in QED for simplicity. Here m refers to the mass of any fermion with electric charge. We are mainly concerned with 
the mass dependence and the differential equation is explicitly given in terms of mass. In the more interesting case of 
symmetry breaking masses, similar analysis should also be feasible and will be studied in the future. 



1 We note that differentiating a Feynman amplitude with respect to external parameters amounts to physically probing the 
quantum process defined by this amplitude, as such operation (d Pli and or 8m) automatically leads to the insertion of elementary 
vertices. Therefore, in this approach a superficially ill-defined amplitude is determined through its response to physical probes 
or 'measurements' plus physical boundary conditions. In conventional renormalization programs, such kind of procedure is 
effected through confronting the intermediately renormalized amplitudes with experiments to obtain physical parametrization 
[10]. 
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This amplitude satisfies the following well-defined inhomogeneous differential equation in any gauge invariant reg- 
ularization (GI) or in the complete underlying theory, 

d m W(p,-p,m) = -ie 2 f (d*k) GI tr{^(—}- )V^+7 M 7-] ^(i^—) 2 }- (2) 

J y /< ~ t1 p + tp-m ft — m p + lfi-m ft — m 

Moreover it satisfies the following well defined equation without regularization at all 
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(d m fW( P ,-p,m) = -^Y,Cf J d*k trm ^^J +W'ij^r 1 } (3) 

where Cf is the combinatorial factor arising from the differentiation operation. 

It suffices to demonstrate our points with Eq. (2). Factorizing out the gauge invariant projector (g^p 2 —p^p") we 
arrive at the following natural equation 

« rr/ 2 \ e 2 f 1 x(l - x)m 

7T 2 J Q m z — x(l - x)p z 

and the solution to this equation reads 

H(p 2 , m; C) = ^ J 1 dxx(l x) In m * ~ x )p* (5) 

with C being the arbitrary integration constant to be fixed somehow. Then we have 



W[p. -p.nr.C) =^(g» v p 2 -p» P nJ o dx (x-x 2 )ln m2 ^ ^ . ((,) 



Here the main conclusion we can draw about C is that it must be gauge invariant and mass independent. 

Now we observe that not all renormalization prescriptions preserve Eq. (4). In all mass independent prescriptions, 
C is fixed so that the renormalized amplitude still satisfies Eq. (4). However, in the MOM schemes, C is defined in 
the following way, 

n(p 2 , m; Cmom) \p>=-v? = Cmom = m 2 + x(l - x)fi 2 , (7) 
which leads to the violation of Eq. (4) , 
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with 5((j?,m) = — dx m -?^^}™^^2 . This is the main point of our concern, as such violation might lead to 
problems in practical applications, especially in the complicated case like SM and its extensions. The spirit followed 
here is the same as is followed in [3,5,6] where physical or novel requirement like gauge invariance are guidelines. In 
these literature the pole mass renormalization for unstable particles is therefore advocated. 

There is a subtle point concerned with the pole mass renormalization for the massless photon, which is just the 
on-shell scheme, as n(p 2 , m; C)| p 2 =0 = — > C = m 2 . Thus it seems that the pole mass definition (or the on-shell 
definition) for photons also violates Eq. (4) . To this end we note that in this case the physical scale or trace anomaly 
[11] is defined in the following way 

e 2 2a 

{p 2 d p2 + m 2 d m 2 - 7A }[1 + n(p 2 , m; C)] = 0, => (3(a) = lA = -2Cd c Ti(p 2 , m; C)\ c=m 2 = — = — (9) 

That is, the mass also play the role of the running scale and the extra term in Eq. (8) is now just the gauge invariant 
trace anomaly with physical significance-the quantum mechanical violation of naive scaling law. Then we see that 
the differential equations could be modified by physical trace anomaly. From now on we consider the violation of the 
differential equations modulo trace anomaly. 

In the more complicated gauge theories such as SM, due to the absorptive (threshold effect) part of the self-energy 
of the unstable particles, mass dependent schemes could lead to gauge dependence and inconsistencey [4], i.e., in such 
theories, the mass dependent schemes were contaminated with severe gauge dependence. We expect that in such 
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theories the novel differential equations proposed here would also be severely violated modulo gauge invariant trace 
anomalies. Thus from the differential equation point of view, combining with the recent studies [2-6] severe gauge 
dependence and other defects, the mass dependent schemes should be avoided in computing SM radiative corrections. 
Theoretically we could not perform a complete sum of the perturbation series to test if such defects automatically 
disappear at all and practically we have to employ 'good' schemes that are free from the above mentioned defects 
before going to experimental data. 

In addition, we should note that the pole mass renormalization for heavy quarks does satisfy the differential 
equations in terms of lagrangian masses modulo trace anomalies, as it is gauge invariant and IR finite [6,12] according 
to the definition m e ff(p 2 , fh)\ p2=M 2 ( = M po i [12]. To see this point we simply note that m e // is a physical 'form 

factor' that is scheme independent [13], thus its dependence upon Lagrangian mass is the same as that in the mass 
independent schemes. This argument also applies to the case of massive bosons, see [2]. 

We can understand the above issues from the underlying theory point of view. Suppose we could compute everything 
from the complete underlying theory, including the 'low energy' amplitudes. Then in the 'low energy' limits, the 'low 
energy' amplitudes must be definite functions (including definite pieces of local terms that arise from the this limit) 
in terms of the 'low energy' physical parameters. There must also be a reference scale for specifying the relevant 'low 
energy' processes. Then the one loop photon polarization should take the following form 

W(p\m u . t .;{a})\ {a} ^ = ^ (g^p 2 - p^p") / dxx{\ - x) In i (10) 

n Jo /Vt. 

with the subscript 'u.t.' referring to the underlying theory and {a} denoting the fundamental constants of the 
underlying theory. It is natural to identify the 'low energy' parameters like m u . t . with the Lagrangian parameters 
of QFT (bare but finite). The scale [i u .t., being process specific and physical and hence renormalization group (RG) 
and renormalization scheme (RS) invariant, is not the running scale though it appears in the same place. Then 
any renormalization prescriptions that are equivalent to such physical parametrization is acceptable. However, the 
difficulty is, we do not know how to identify such process dependent scale theoretically before the complete underlying 
theory is known, therefore we must resort to other means like experimental data and some reasonable or physical 
requirements. Before the complete theory is found, we could not exclude the possibility that there might be schemes 
that are physically inequivalent to the complete theory definitions at all. 

Moreover, from the differential equation point of view, the scheme invariance of certain quantities in QFT originally 
established within the mass independent schemes [1] could be promoted to a full 'invariance' for all the schemes that 
satisfy the differential equations in terms of all physical parameters. Otherwise, if the mass dependent schemes also 
qualified for applications, the significance of such 'invariance' established within the mass independent schemes is 
diminished in the cases where masses could not be neglected. In a sense we put forward an argument for advocating 
mass independent schemes, and this point of view could draw supports from the recent literature on renormalization 
schemes, though the relation between our approach here and those in the literature is unclear at this moment. 

Before closing the presentation we would like to discuss the issue of heavy particles decoupling where the MOM 
and other mass dependent schemes are thought to be advantageous over the mass independent schemes due to its 
good decoupling behavior. We first note that this advantage is only in the context of RGE. Then we also remind that, 
as we will demonstrate below, in the context of Callan-Symanzik equations [14] (CSE), which describe the full scale 
behavior, the decoupling of heavy fields [15] is achieved in the same way in all schemes. 

Again we illustrate it with a simple model, QED with a massive fermion in addition to ni massless fermions. In a 
mass independent scheme the Callan-Symanzik equation reads, 

{\d x - (iad a + 7r - D r }T((Xp),m, a, = -iT e {{\p), to, a, M ) (11) 

where = [1 + r y m ]m r ip^p, 0, -fr and -f m are mass independent functions of the renormalized coupling a and all 
quantities are renormalized ones. At lowest order, = |^[n; + 1]. 

When the mass goes to infinity, 0, 7 on the left hand side of CSE does not change at all, but in the meantime, the 
inhomogeneous term of CSE does not vanish, 

- J F e ((p),m,a, A i)|„ woo - (A/3a9 Q - A 7r )r ((p),a,M), (12) 
=}► {Xd x - a ad a + 7r ;0 - D r }T ((Xp), a, M ) = 0. (13) 

Here 0o = + A/3, 7r ; o = 7r + A 7 r with the delta contributions coming from the mass insertion part in the infinite 
mass limit which will cancel the heavy field's contribution to and 7. The subscript means that the heavy particle 
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is removed. The generalization to other theories with boson masses is an easy exercise. From Eq. (13) we see that 
the decoupling of heavy particles is realized in a natural way in the contexts of Callan-Symanzik equations. 

To verify the above deduction it is enough to demonstrate Eq. (12) at the lowest order which is closely related to 
the observation that heavy particle limit provides a convenient algorithm for calculating trace anomalies [16]: 

-i^^J" J")| m ^oo - ^(/r-?V) => m(l + 7mWHm^oo = \ ApF^F^, (14) 

with = —ie'ijj'j^ip and A/3 = — |f = Aj^. When translated into Callan-Symanzik equations Eq. (14) is just 
Eq. (12). The cancellation of the heavy particle contributions is obvious since (3 + A/3 = |f [n/ + 1] — |f = §f«;- 
While in the MOM like schemes the Callan-Symanzik equation reads, 

{Xd\ - pMOMOida + Jc-mom - D T }T mom({Xp), m, a, /z) = -iVf {OM ({\p), m, a, fj,) (15) 

with the beta function etc. being defined as Pmom = [pd^ + m(l + r ) m ;MOM)d m ]oi, • ■ • , in contrast to the RGE 
definition: /?mom = A^a, • • • , due to the mass dependence of the renormalization constants 2 . The resulting (3mom 
also exhibits non-decoupling feature as in the mass independent schemes. For example, at the lowest order, from the 
definition given above, a heavy particle's contribution to (3 at lowest order is mass independent 

i8mom ; m = (K + m(l + 7 m ;MOM)9m)^-2 / dxx(l - x) In Cmom = (16) 

Without any doubt such 'non-decoupling' term is cancelled by the contribution from the decoupling limit of the 
inhomogeneous term in CSE, that is, just like in Eq. (12), we have 

_ *r MOM ((p),m, a,[l,)\ m ->oo = {ApMOMO-da - Ajr:MOM)T A4OM:0((p) , a, fl) , (17) 

with A[3mom = — §f j which exactly cancels the heavy particle's contribution to the (3, 7 etc. Thus we proved that 
Eq. (12) is also true in the MOM schemes at the lowest order. It is well known that the first loop order beta of RGE 
in mass independent schemes differs from that in the MOM schemes [13]. While in the context of Callan-Symanzik 
equation the beta function is the same in all schemes at one loop order (C.f. Eq. (16)), which in turn implies that 
the same decoupling mechanism works in both mass independent and mass dependent schemes in the context of 
Callan-Symanzik equation. Of course, in the preceding derivation it is assumed that there is no such trouble as gauge 
dependence and/or any form of inconsistency for the mass dependent schemes, i.e., we consider the theories that are 
not beset with the troubles that afflicts the electroweak sector of SM. 

In summary, we revisited the issue of renormalization schemes that is intensely discussed in the literature of Standard 
Model. The analysis, basing on a set of natural differential equations, seems to favor the mass independent schemes 
against the mass dependent ones, supporting the point of view of the recent literature concerning the renormalization 
of CKM matrix and the unstable sectors of SM. It is also demonstrated that in the context of Callan-Symanzik 
equation all the schemes facilitate the decoupling of heavy particles in the same way. Further investigation in this 
direction seems worthwhile. 



ACKNOWLEDGEMENT 



The author is grateful to W. Zhu for useful discussions. This work is supported by the National Nature Science 
Foundation of China under Grant No. 10075020 and No. 10502004. 



2 In the standard inhomogeneous form of Callan-Symanzik equation, the mass operators inserted vertex functions appear in the 
other side of the equation and the definitions of f3, 7 in the MOM-like mass dependent schemes must include m(l + ~/m;MOM)d m 
to account for the 'full running' of the renormalization constants. But in the alternative homogeneous form of Callan-Symanzik 
equation (CSE), i.e., {\d\ — /3ad a + m(l + j m )dm + 7r — ^r}r((Ap), m, a, n) = 0, the definitions of the /?, 7 etc. are the same 
as in RGE. 
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